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One believes there is huge amount of Dark Matter particles in our Galaxy which 
manifest themselves only gravitationally. There is a big challenge to prove their 
existence in a laboratory experiment. To this end it is not sufficient to fight only 
for the best exclusion curve, one has to see an annual recoil spectrum modulation 
— the only available positive direct dark matter detection signature. A necessity to 
measure the recoil spectra is stressed. 

PACS: 95.30.-k, 95.35.+d, 14.80.Ly, 12.60.Jv 
Galactic Dark Matter (DM) particles do not emit (or reflect) any detectable electromagnetic 

O-i 

radiation and manifest themselves only gravitationally by affecting other astrophysical objects. 



According to the estimates based on a detailed model of our Galaxy [l] the local density 
of DM (nearby the solar system) amounts to about pj^jj — 0.3 GeV/cm 3 ~ 5 ■ 10~ 25 g/cm 3 
(see also recent reviews 0, [j]). The local flux of DM particles x is expected to be — 

— — • 10 5 cm~ 2 s _1 , where m y is the DM particle mass. This value is often considered as 

m x 

a promising basis for direct laboratory dark matter search experiments. 

The problem of the DM in the Universe is a challenge for modern physics and experimental 
technology. To solve the problem, i.e. at least to detect the DM particles, one simultaneously 
needs to apply the front-end knowledge of modern Particle Physics, Astrophysics, Cosmology 
and Nuclear Physics and to develop and use over long time extremely high-sensitive experi- 
mental setups and complex data analysis methods (see, for example, recent discussion in |4|). 

Weakly Interacting Massive Particles (WIMPs) are among the most popular candidates for 
the relic DM. These particles are non-baryonic and there is no room for them in the Standard 
Model of particle physics (SM). The lightest supersymmetric (SUSY) particle (LSP), neutralino 
(being massive, neutral and stable), is currently often assumed to be a favorite WIMP dark 
matter particle. 

The nuclear recoil energy due to elastic WIMP-nucleus scattering is the main quantity to be 
measured by a terrestrial detector in direct DM detection laboratory experiments 5j. Detection 
of the very rare events of such WIMP interactions is a quite complicated task because of very 
weak WIMP coupling with ordinary matter. The rates expected in the SUSY models range 



from 10 to 10~ 7 events per kilogram detector material a day {6-[l3|. Moreover, for WIMP masses 
between a few GeV/c 2 and 1 TeV/c 2 , the energy deposited by the recoil nucleus is less than 
100 keV. Therefore, in order to be able to detect a WIMP, an experiment with a low-energy 
threshold and an extremely low radioactive background is required. Furthermore, to certainly 
detect a WIMP one has to unambiguously register some positive signature of WIMP-nucleus 
interactions (directional recoil or annual signal modulation) ?} Q]- This means one has to 
perform a stable measurement with a detector of large target mass during 3-5 years under 
extremely low radioactive background conditions. There are also some other complications 



discussed recently in 



4]- 
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16] has certainly observed the 



Till now only the DAMA (DArk MAtter) collaboration (4, 
first evidence for the DM signal due to model-independent registration of the predicted annual 
modulation of specific shape and amplitude due to the combined motions of the Earth and the 
Sun around the galactic center [14J. This experiment has released a total exposure of 1.17 t 
x yr over 13 annual cycles, obtaining positive model-independent evidence for the presence of 
DM particles in the galactic halo at 8.9 a CL. 



Although there are other experiments like EDELWEISS, CDMS, XENON, CRESST, etc, 
which give sensitive exclusion curves, no one of them at present has the sensitivity to look for 
the modulation effect. Due to the relatively small target masses and short running times these 
experiments are unable to see a positive annual modulation signature of the WIMP interactions. 
Unfortunately, some other experiments with targets of much larger mass (mostly Nal) were also 



unable to register the positive signature due to not good enough background conditions |17H19|. 

Despite the strong and reliable belief of the DAMA collaboration in the observation the 
annual modulation signature, it is obvious that such a serious claim should be verified by at 
least another one completely independent experiment. 

If one wants to confirm (more important, if one wants to reject) the DAMA result, one 
should perform a new experiment which would have the same or better sensitivity to the annual 
modulation signature (and also it would be reasonable to locate this new setup in another low- 
background underground laboratory). In particular, search for the modulation could be carried 
out by new-generation experiments with high purity germanium detectors of large enough 
mass, perhaps, both with spin 73 Ge and spinless natural Ge 20|. It is interesting that recently 



the CoGeNT experiment with a germanium detector has reported some preliminary positive 
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indication of the annual modulation 



2l|. 



Together with necessary figthing against backgrounds, the main direction in development 
of new-generation DM detectors concerns remarkable enlargement of the target mass to allow 
observing these positive signatures and thus detecting DM and proving or disproving the DAMA 



claim. In particular, an enlarged version of the EDELW 



detectors 



22] tog ether with, perhaps, SuperCDMS 
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USS setup with 40 kg bolometric Ge 



24J, as well as the enlarged ZEPLIN 



25] or KIMS 26] experiments might become sensitive to the annual modulation in the future. 



To estimate the expected direct detection rate for these WIMPs (in particular, neutralinos) 
any SUSY-like model or some measured data, for example, from the DAMA experiment 271 ]. 
can be used. On this basis the WIMP-proton and WIMP-neutron spin cr$£(0) and scalar 
agj n (0) cross sections at zero-momentum transfer can be calculated (see the Appendix). These 
calculations are usually compared with measurements, which (with the only exception of the 
DAMA result) are presented in the form of exclusion curves — upper limits of the cross section 
as functions of the WIMP mass. In the case of non-observation of any DM signal the exclusion 
curve simply reflects the sensitivity of a given direct DM search experiment and potentially 
allows one to constrain some version of the SUSY-like theory, if the curve is sensitive enough. 
Therefore the best exclusion curve is currently a clear aim of almost all dark matter search 
experiments (DAMA/LIBRA and CoGeNT are the only exceptions). The main competition 
between the experiments is in the field of these exclusion curves. 

Before 2000 all exclusion curves were evaluated mainly in the one-coupling dominance ap- 
proach (when only one cross section limit was defined from measurements for fixed WIMP 
mass), which gave slightly pessimistic (for spin-non-zero target experiments), but universal 
limits for all experiments. One would say that the competition between the DM experiments 
was honest. The predictions from SUSY-like models were in general far from being reached by 
the data. 



Mainly after the paper [28] was published in 2000 (and as well after the DAMA evidence 
15] ) a new kind of exclusion curves appeared. In particular, for the first time these curves were 
obtained for the spin- dependent WIMP-nucleon cross section limits when non-zero subdomi- 
nant spin WIMP-nucleon contributions were also taken into account (29), 0]- This procedure 
obviously improved the quality of the exclusion curves. Therefore a direct comparison of the old- 
fashioned exclusion curve with the new one could in principle bring one to a wrong conclusion 
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about better sensitivity of more recent experiments. There is generally possible incorrectness 
in the direct comparison of the exclusion curves for the WIMP-proton(neutron) spin-dependent 
cross section obtained with and without the non-zero WIMP-neutron (proton) spin-dependent 
contribution. Furthermore, the above-mentioned incorrectness concerns to a great extent the 
direct comparison of the spin-dependent exclusion curves obtained with and without non-zero 



spin-independent contributions 



15 



311 ] . Taking into account both spin couplings a p and 



but ignoring the scalar coupling cq (see the Appendix for definitions), one can easily arrive at 
a misleading conclusion especially for not very light target nuclei when it is not obvious that 
(both) spin couplings dominate over the scalar one. To be consistent, one has to use the mixed 



spin-scalar coupling approach as was first proposed by the DAMA collaboration [15|, |31|, 132 ] . 

This approach was used in 33| to demonstrate, by the example of the HDMS experiments 
with natural Ge and with the neutron-odd group high-spin isotope 73 Ge {34] . 35], how one can 
strongly improve the exclusion curves. The approach allowed both upper limits for the spin- 
dependent <Tg£^ and spin-independent a sl cross sections of the WIMP-nucleon interaction to 



be simultaneously determined from the experimental data. In this way visib 



magnitude) improvement in the form of the exclusion curves was achieved 33] relative to the 



e (one order of 



1 



traditional one-coupling dominance scheme used previously for the same setup 

As a by-product of the approach, there are correlations (first mentioned in [37]) between 
the mesured upper limits <Tg D and a sl , which can be considered as a new requirement - 
for any fixed WIMP mass m x one should have <7 SI (theor.) < cr gi (exp.) and <jg D (theor.) < 
<7g D (exp.) simultaneously, provided that crg^g^(theor.) are calculated in some underlying SUSY- 
like theory. 

It is important to note that without proper knowledge of the nuclear and nucleon structure 
it is not possible to extract reliable and useful information (at least in the form of these <7g D and 
crgj cross sections) from direct DM search experiments. However, astrophysical uncertainties, 
in particular the DM distribution in the vicinity of the Earth 38l443j , make it far more dificult 
to interpret the results of the DM search experiments. At the moment, to have a chance to 
compare sensitivities of different experiments, people adopted a common truncated Maxwellian 
DM particle distribution, but nobody can prove its correctness. In the case of undoubted direct 
DM detection one can make some conclusions about the real DM particle distribution in the 
vicinity of the Earth. 
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Furthermore, almost by definition (from the very beginning), a modern experiment aiming 
at the best exclusion curve is doomed to non-observation of the DM signal. This is due to the 
fact that a typical expected DM signal spectrum exponentially drops with recoil energy and it 
is practically impossible to single it out from the background non-WIMP spectrum of a typical 
(semiconductor) detector. 

In fact, one needs a clear, so-called "positive" signature of interactions between WIMP 
particles and target nuclei. Only exclusion curves are not enough. Ideally, this signature 
should be a unique feature of such an interaction [44 1. 

There are some typical characteristics of WIMP particle interactions with a nuclear target 
which can potentially play the role of these positive WIMP signatures 45j. First of all, WIMPs 



produce nuclear recoils, whereas most radioactive backgrounds produce electron recoils. Never- 
theless, for example, neutrons (and any other heavy neutral particle) can also produce nuclear 
recoils. There are also proposals which rely on WIMP detection via electron recoils 46l. 147 ] . 



Due to the extremely rare event rate of the WIMP-nucleus interactions (the mean free path 
of a WIMP in matter is of the order of a light year), one can expect two features. One is 
that the probability of two consecutive interactions in a single detector or two closely located 
detectors is completely negligible. Multiple interactions of photons, gamma rays or neutrons 
under the same conditions are much more common. Therefore only non-multiple interaction 
events can claim to be from WIMPs. The other feature is a uniform distribution of the WIMP- 
induced events throughout a detector. This feature can also be used in the future to identify 
background events (from photons, neutrons, beta and alpha particles) in rather large-volume 
position-sensitive detectors. 

The shape of the WIMP-induced recoil energy spectrum can be predicted rather accurately 
(for given WIMP mass, fixed nuclear structure functions, and astrophysical parameters). The 
observed energy spectrum, claiming to be from WIMPs, must be consistent with the expecta- 
tion. However, this shape is exponential, right as it is the case for many background sources. 

Unfortunately, the nuclear-recoil feature, the non-multiple interaction, the uniform event 
distribution throughout a detector, and the shape of the recoil energy spectrum could not be 
the clear "positive signature" of the WIMP interactions. It is believed that the following three 
features of WIMP-nucleus interaction can serve as a clear "positive signature" . 

The currently most promising, technically feasible and already used (by the DAMA collabo- 
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ration) "positive signature" is the annual modulation signature (see the Appendix). The WIMP 
flux and its average kinetic energy vary annually due to the combined motions of the Earth and 
the Sun relative to the galactic center. The impact WIMP energy increases (decreases) when 
the Earth velocity is added to (subtracted from) the velocity of the Sun. The amplitude of the 
annual modulation depends on many factors — details of the halo model, mass of the WIMP, 
the year-averaged rate (or total WIMP-nucleus cross sections), etc. In general, the expected 
modulation amplitude is rather small fly , ll5|, |3l| and to observe it, one needs huge (at best 
tonne scale) detectors which can continuously operate for 5-7 years. Of course, to reliably use 
this signature one should prove the absence of annually modulated backgrounds. 

Another potentially promising positive WIMP signature is connected with the possibility 
of measuring the direction of the recoil nuclei induced by a WIMP. In these directional recoil 
experiments it is planned to measure the correlation of the event rate with the Sun's motion 
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49]. Unfortunately, the task is extremely complicated 5(M54 1 . 



The third well-known potentially useful positive WIMP signature is connected with the 
coherence of the WIMP-nucleus spin-independent interaction. Due to a rather low momentum 
transfer, a WIMP coherently scatters by the whole target nucleus and the elastic cross section 
of this interaction should be proportional to A 2 , where A is the atomic number of the target 
nucleus. Contrary to the A 2 behavior, the cross section of neutron scattering by nuclei (due 
to the strong nature of this interaction) is proportional to the geometrical cross-section of the 
target nucleus (A 2 ^ 3 dependence). To reliably use this A 2 signature, one has to satisfy at least 
two conditions. First, one should be sure that the spin-independent WIMP-nucleus interaction 



real 



33 



y 



dominates over the relevant spin-dependent interaction. This is far from being obvious 
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57j . Second, one should, at least for two targets with a different atomic number A, rather 
accurately measure the recoil spectra (in the worst case integrated event rates) under the same 
background conditions. Currently, this goal looks far from being achievable. Developing further 
the idea of this third signature, one can also consider as a possible extra WIMP signature an 
observation of the similarity (or coherent behavior) of measured spectra at different (also non- 
zero spin) nuclear targets. This possibility relies on rather accurate spin structure functions for 
the experimentally interesting nuclei 58L l59|. 

Ideally, in order to be convincing, an eventual DM signal should combine more than one of 



these positive DM signatures 



44 



45]. 
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In the case of currently very promising event-by-event active background reduction tech- 
niques (like in the CDMS, EDELWEISS and XENON experiments), one inevitably needs clear 
positive WIMP signature(s). Without these signatures one can hardly convince anybody that 
the final spectrum is saturated by WIMPs. Furthermore, with the help of these extra signatures 
and on the basis of measured recoil spectra one can estimate the WIMP mass 60|, l61 |. 



It is known (see, for example, discussions in 62|, |63|) that a proof of the observation of a DM 



signal is an extremely complicated problem. As pointed out above, on this way an interpretation 
of measurements in the form of exclusion curves helps almost nothing. Of course, an exclusion 
curve is at least something from nothing observed. It allows a sensitivity comparison of different 
experiments and therefore allows deciding who at the moment is the best 'excluder'. But, for 
example, supersymmetric theory is, in general, very flexible, it has a lot of parameters, and 
one hardly believes that an exclusion curve can ever impose any decisive constraint on it. 
Furthremore, almost all experimental groups presenting their exclusion curves try to compare 
them with some SUSY predictions. It is clear from this comparison that there are some domains 
of the SUSY parameter space, which are now already excluded by these exclusion curves. What 
is remarkable, however, is that nobody yet has seriously considered these constraints for SUSY. 
The situation is much worse due to the already mentioned famous nuclear and astrophysical 



uncertainties involved in the evaluation of the exclusion curves 
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71| . This is why it does 



not look very decisive (or wise) to use very refined data and methods (nuclear, astrophysical, 



numerical, statistical 



72| , etc) and spend big resources fighting only for the best exclusion 



curve. This fighting could only be accepted when one tries to strongly improve the sensitivity 
of a small detector with a view of using many copies of it in a huge detector array with a total 



tonne-scale mass 



20]. 



There are remarks concerning comparison of results from DM search experiments with pas- 
sive (off-line) background reduction (like DAMA) and from experiments with active (on-line) 
background reduction (like CDMS, XENON, ZEPLIN, etc). First, it was demostrated Q- 
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33] that any extra positively defined background-like contribution to the spectra improves 
the extracted (upper limit) values of the cross section. Next, within the passive background 
reduction scheme the measured spectrum is not affected by hardware or software influence dur- 
ing the data taking. Further background reduction can be done off-line on the basis of careful 
investigation of the spectrum itself or, for example, with the help of the pulse shape analysis. 



s 



In this case the extracted background contribution is under control and well defined. On the 
other side, within the active background reduction approach the measured spectrum already 
contains results of this active reduction influence on the data taking process. In this case it 
is not simple to hold under control the real level of extracted on-line background contribution 
which can easily be overestimated (see, for example, the relevant discussion in |4j). Therefore, 
due to this obvious difference, a direct comparison of exclusion curves from experiments with 
passive and active background reductions could be, in principle, rather misleading. 

Finally, it seems that at the level of our present knowledge the DM problem could not 
be solved independently of other related problems (proof of SUSY, astrophysical dark matter 
properties, etc). Furthermore, due to the huge complexity of the DM search (technical, physical, 
astrophysical, necessity for positive signatures, etc), one should deal with the DM problem 
boldly using a reliable model-dependent framework — for example the framework of SUSY, 
where the same LSP neutralino should be seen coherently or lead to effects in all available 
experiments (direct and indirect DM searches, rare decays, high-energy searches at LHC, etc). 
Only if such a SUSY framework leads to a specific and decisive positive WIMP signature, this 
could mean a proof of SUSY and simultaneous solution of the dark matter problem. In some 
sense, this SUSY framework can serve as a specific and very decisive positive WIMP signature. 



SUMMARY 



A physical reason to improve an exclusion curve is usually an attempt to constrain a SUSY- 
like model. Unfortunately, this is almost hopeless due to the huge flexibility of these models 
and the inevitable necessity of having extra information from other SUSY-sensitive observables 
(for example, from LHC). At the present and foreseeable level of experimental accuracy, sim- 
ple fighting for the best exclusion curve is almost useless either for real DM detection or for 
substantial restrictions for SUSY. 

One should inevitably go beyond an exclusion curve. New generations of DM experiments 
right from their beginning should aim at detection of the DM particles. This will require 
development of new setups, which will be able to register positive signatures of the DM particles 
interactions with nuclear targets. 

One should try to obtain a reliable recoil energy spectrum. First, very accurate off-line 
investigation of the measured spectrum allows one to single out different non-WIMP background 
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sources and to perform controllable background subtractions. Second, the spectrum allows one 
to look for the annual modulation effect, the only currently available positive signature of DM 
particle interactions with terrestrial nuclei. This effect is not simply a possibility (among many 
others) of rejecting background (as claimed again receintly in |3(), but it is a unique signature 
which reflects the inner physical properties of the DM interaction with matter. It is a very 
decisive and eagerly welcomed feature, which is inevitable for the laboratory proof of the DM 
exsitence. 

This letter was written in connection with Prof.D.I.Kazakov's 60th birthday and contains 
updated key messages from the extended review "Direct Search for Dark Matter — Striking 
the Balance — and the Future" (20 ]. 

Preparation of this work was supported by the grant of the Ministry of Education and 
Science of the Russian Federation (contract 12.741.12.0150). 



APPENDIX 



The nuclear recoil energy E R is measured by a proper detector deep underground (Fig. [p. 
The differential event rate in respect to the recoil energy (the spectrum) is the subject of 
the measurements. The recoil spectrum produced from WIMP-nucleus scattering in a target 
detector is expected to show the annual modulation effect due to the Earth's motion around the 
Sun 14]. The velocity of the Earth relative to the Galaxy is Vs{t) = v$ +vo cos 7 cos u(t — t ), 
where vs is the Sun's velocity relative to the Galaxy (ps = 232 km/s), vo is the Earth's orbital 
velocity around the Sun (vq = 30 km/s) and 7 is the angle of inclination of the plane of 
the Earth's orbit relative to the galactic plane (7 = 60°). One has u = 2n/T (T = 1 year) 
and the maximum velocity occurs at day to — 155.2 (June 2). The change in the Earth's 
velocity relative to the incident WIMPs leads to a yearly modulation of the scattering event 
rates of about 7%. It is convenient to introduce a dimensionless variable 77 = ve/vq, then 
T)(t) = 7/0 + A 7/ cos 00 (t — to) where the amplitude of the modulated part (Arj ~ 0.07) is small 
compared to the annual average r] ~ 1.05. Within this framework, the expected count rate of 
WIMP interactions can be written as 



81 

S[ V (t)\ ~ % ] + f- 



Ar] cos u(t — to) = Sq + S m cos w(t — to) 

no 
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FIG. 1: Detection of dark matter (WIMPs) by elastic scattering from target nuclei in the detector. 
Due to the expected annual modulation signature of the event rate J2]), the Sun-Earth system is a 
particularly proper setup for successful direct DM detection. From [201 ] . 



where So is the constant part and S m is the amplitude of the modulated signal. Both parts 
of the event rate S and S m depend on the target nucleus (A,Z), WIMP (or neutralino x) 
mass m x , density p^^j, velocity distribution of the WIMPs in the solar vicinity f(v), and cross 
section of WIMP-nucleus elastic scattering (see, for example, 
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The differential event rate per unit mass of the target material has the form 



S(t) = 



dR 



n fmai d(T"^ 



dq 2 



(2) 



Assuming that WIMPs are the dominant component of the DM halo of our Galaxy, one has p x = 
Piocli- The nuclear recoil energy E R = q 2 /{2MA) is typically about 10 _6 m x , N T is the number 



density of target nuclei with mass Ma, v n 



1 /2 

600 km/s, and v min = (M A E R /2p 2 A ) 1 is 



the minimal WIMP velocity which still can produce the recoil energy E^. The WIMP-nucleus 



differential elastic scattering cross section for spin-non-zero (J 7^ 0) nuc 



(spin- independent, or SI) and axial (spin-dependent, or SD) terms [75|, [76] 



ei contains coherent 



da' 
~dq~ 2 



[v,q 



v 2 (2J + l) v 2 (2J+l 



4d(0) 



4fi 2 A v 2 * SD 



4(0) p2 ,^ 



4p> 



(3) 



The normalized (Fg D gi(0) = 1) finite- momentum-transfer nuclear form-factors Fg DSl (q 2 
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*S'sD,Sl('? 2 ) 
^SD,Sl(0) 



can be expressed in terms of the nuclear structure functions as follows 75|, |76| : 



S&(q) = E |(J||C L (g)||J)r^|(J||C (g)||J)| 2 , 

L even 

StM = E (K^||r^ 5 (g)||iV)| 2 + |(iV||£|(g)||iV)| 2 ). 



(4) 



Lodd 



The explicit form of the transverse electric T el5 (q) and longitudinal C 5 (q) multipole pro 



of the axial vector current operator and the scalar function Ci(q) can be found in 58, 
For q = the nuclear SD and SI cross sections can be represented as 



59 



ections 



75 



76|. 



4(0) 



4/4S SI (0) ii\ 



f4 



SD 



(2J+1) 

crZ (0) { (S*) cos 9+(S*} sin^} 2 . 



(2 J 

/£ 3 J 



(5) 
(6) 
(7) 



Following Bernabei et al. 15|, |37|], the effective spin WIMP-nucleon cross section cxgp(O) and 
the coupling mixing angle 9 were introduced, 



U 4 

*§£(o) = ^ 

7T O 



a 2 + a 2 n 



tan# = — ; 



a 



SD 



pn 2/i n pn • 2 n 

ag D ■ cos 9, a SD = crg D • sin 9. 



(9) 



Here, /x^ 



m x M A 
m x + Ma 



is the reduced mass of the neutralino and the nucleus, and it is assumed 



that /i 2 = fip. The dependence on effective WIMP-quark (in SUSY neutralino-quark) couplings 
C g and A q in the underlying theory 



£es = E (A • • qi^i^q + c q -xx-qq) +■■■ 



(10) 



and on the spin (A|f' n )) and the mass or scalar (f^ m /^ n - ) ) structure of the proton and neutron 
enter into these formulas via the zero-momentum-transfer WIMP-proton and WIMP-neutron 
SI and SD cross sections 



<(0)=4^c 2 , 

71 



^sb (0) = 12— a 

7T 



r — ^ n — vr f(p."). 

c ~~ C — ^QJq ' 

ctp = ^4 ? A g , a n = .Aq ^ . 

9 



(11) 

(12) 
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The factors A q p ' n \ which parameterize the quark spin content of the nucleon, are defined as 
2A^"' P V' = (p, s\ip q 'y fJ ''j5ip q \p, s)(p,n)- The quantity (Spix) denotes the total spin of protons 
(neutrons) averaged over all A nucleons of the nucleus (A, Z) 



Kn)) = (A\S* n) \A) = (A\J2si (n) \A). 



(13) 



The mean velocity (v) of the relic neutralinos of our Galaxy is about 300 km/s = 10~ 3 c. 
Assuming q ma , x R <C 1, where R is the nuclear radius and g max = is the maximum of the 

momentum transfer in the process of the xA scatterin g, t he spin-dependent matrix element 
takes a simple form (zero momentum transfer limit) 77J, |78| 



M = C(A\a p S p + a n S n \A) ■ s x = CA(A\J\A) ■ s x . 



Here, s x denotes the spin of the neutralino, and 



A 



(N\a p S p + a n S n \N) _ (N\(a p S p + a n S n ) ■ 3\N) _ a p {S p ) a n (S n ) 



(14) 



(15) 



(iV|J|A0 J(J+1) J J 

The normalization factor C involves the coupling constants, the masses of the exchanged bosons, 
and the mixing parameters relevant to the LSP, i.e., it is not related to the associated nuclear 
matrix elements (tq| . In the limit of zero momentum transfer q = the spin structure function 
in (j3J) reduces to the form 

S A (0) = ±-\(A\\ £ i(a + air>,||A)| 2 = ^±^J(J + 1)A 2 . 
For the most interesting isotopes either (S^ ) or (S^) dominates ((Sf^,) <C (S^)). 



The differential event rate ([2]) can be also given also in the form 
dR(E R ) 
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55] 



dE R 
K sl (E R ,m x ) 

K SD (E R ,m x ) 

Bsi,sd(E r ) 



k si (E r , m x ) (T S i + k S d(£r, m x) °SD- 
N T ^B sl (E R ) 

Nt^t — ^-Ssd(-E'r) 



(16) 



2m x [L 2 p 



4 J + l 
3 J 



((S p ) cos 9 + (S n ) sin i 



(17) 



<^ 2 > 



The dimensionless integral I(E R ) is a dark-matter-particle velocity distribution correction 



I(E ¥ 



{ft f m 

(v) 



dx 



3 + 2r/ 2 



2 y/Tr(l + 2r] 2 )eri(r]) + 2r)e~^ 



■[erf(x min + rf) - erf(x min - rf)], 



(18) 
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where WIMPs in the rest frame of our Galaxy are assumed to have a Maxwell-Boltzmann 
velocity distribution, the dimensionless Earth speed with respect to the halo r], is used, and 



x 



3M A E_ 



The error function is erf(x) = — = / dte . The velocity variable is 

7T JO 



4 \i\v 2 

the dispersion v ~ 270km/c. The mean WIMP velocity (v) = \f\v. Integrating the differential 
rate fl2]) from the recoil energy threshold e to some maximal energy e, one obtains the total 
detection rate R(e, e) as a sum of the SD and SI terms 

R(e,e) = R S i(e, e) + R SD (e, e) = J dE R n sl (E Rj m x ) a S i + J dE R K SD (E R ,m x ) a SB . (19) 

To accurately estimate the event rate R(e,e), one needs to know a number of quite uncertain 
astrophysical and nuclear structure parameters as well as the very specific characteristics of the 



qR becomes non-negligible and the finite momentum transfer 



experimental setup 

As m x increases, the product 



limit must be considered 



58, 
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75 



76 



781 ] . With the isoscalar spin coupling constant a 



0"n + o, p and the isovector spin coupling constant a\ = a p — a n , one can split the nuclear 
structure function S A (q) into a pure isoscalar term, S^q), a pure isovector term, S^q), and 
an interference term, S^q), in the following way: 



S A (q) = a 2 S A (q) + afS^q) + ao^S^q). 



2 a A 



(20) 



The relations ,S A (0) = C{J){{S P ) + (S n )) 2 , 5^(0) = C(J)((S P ) - (S n )) 2 , and S£(0) = 

2C(J)((S 2 ) — (S 2 )) with C(J) = — ^— connect the nuclear spin structure function 

S A (q = 0) with the proton (S p ) and neutron (S„) spin contributions averaged over the nucleus 



58]- 



To analysis modern data in the finite momentum transfer approximation it seems reasonable 
to use the formulas for the differential event rate (j2J) as schematically given below 
dR(e, e) 



dE R 
Af(e,e,E R ,m x ) 

r) SI (E R ,m x ) 



W(e, e, E R , m x ) T] sl (E R , m x ) a p sl + r]' SD (E R , m x , u) a 



(21) 



A^ 7 



cp x M A 
2m x n 2 v 



(?max) C 



)I(E R )9(E R - e)6{e - E R ) 



{^s 2 i(^r)}; 



i-4 



2J 



■J (Soo(q) + to 2 Su(q) + u S 01 (q)) > . 



Here the isovector-to-isoscalar nucleon coupling ratio is u = ai/oo- The detector threshold 
recoil energy e and the maximal available recoil energy e (e < E R < e) have been introduced 
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in f lT^j) . In practice, with an ionization or scintillation signal, one has to take into account the 
quenching of the recoil energy, when the visible recoil energy is smaller than the real recoil 
energy transmitted by the WIMP to the target nucleus. 

Formulas (|2ip allow experimental recoil spectra to be directly described in terms of only 
three [80[ (it is rather reasonable to assume ergi(O) ~ crgi(O)) independent parameters (crg I; 0% 
and u) for any fixed WIMP mass m x and any neutralino composition. Comparing this formula 
with the observed recoil spectra for different targets (Ge, Xe, F, Nal, etc) one can directly and 
simultaneously restrict both isoscalar cq (via Ugj) and isovector neutralino- nucleon effective 
couplings ao x- These constraints, based on the nuclear spin structure functions for finite q, will 
impose the most model-independent and most accurate restrictions on any SUSY parameter 
space. Contrary to some other possibilities (see, for example, |15|, [28[), this procedure is direct 
and uses as much as possible the results of the accurate nuclear spin structure calculations. 

It is seen from (jH]) and (liZTj) that the SD cross sections <Tg D and <7g D (or equivalently 
and u = ai/ao) are the only two WIMP-nucleon spin variables which can be constrained (or 
extracted) from DM measurements. Therefore, there is no sense in extracting effective WIMP- 
nucleon couplings a and a n from the data (with "artificial" twofold ambiguity). 

Finally, to estimate the expected direct DM detection rates (with formulas (T5]), ( 1191) or (ED)) 
one should calculate the cross sections <r SI and cr SD (or WIMP-nucleon couplings Co and a p n ) 
within a SUSY-based model or take them from experimental data (if it is possible). 
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